Negative index of refraction in a four-level system with magnetoelectric cross coupling 

and local field corrections 
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This research focuses on a coherently driven four-level atomic medium with the aim of inducing a 
negative index of refraction while taking into consideration local field corrections as well as magneto- 
electric cross coupling, i.e. chirality, within the material's response functions. Two control fields are 
used to render the medium transparent for a probe field which simultaenously couples to an electric 
and a magnetic dipole transition, thus allowing one to test the permittivity and permeability of the 
material at the same time. Numerical simulations show that a negative index of refraction with low 
absorption can be obtained for a range of probe detunings while depending on number density and 
the ratio between the intensities of the control fields. 
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I. INTRODUCTION 

The search for a negative refractive index within the 
visible range of frequencies has attracted a considerable 
amount of interest due in part to the advances made in 
the production of optical metamaterials (see for instance 
[l| and references therein), where the optical properties 
of the medium are determined mainly by its structure 
rather than by its intrinsic properties. Materials with 
a negative index of refraction, also known as left-handed 
materials, have a negative phase velocity with power flow 
in the opposing direction, and have been theoretically 
proven to exhibit phenomena such as subwavelength reso- 
lution, a reverse Doppler effect, and Cherenkov radiation 
to name a few [3, y] . Although experiments have veri- 
fied the existence of negative index materials and a wide 
range of applications [jig, losses remain a major limi- 
tation in the optical regime [3-[y]- Alternatively, it was 
shown recently that negative refractive index can be clec- 
tromagnetically induced in vapors with low absorption 
rates |10l . [ll| thus eliminating the need of manufactur- 
ing the medium. In addition, recent publications [12|-[l8j 
have pointed out that if one considers magnetoelectric 
cross coupling and/or local field corrections, a negative 
index of refraction can be achieved without requiring a 
negative permeability, which so far has been one of the 
largest impediments to overcome in the field. 

Using density matrix theory in order to treat the in- 
teraction of light with matter semiclassically, this pa- 
per studies a dense medium of 4-level atoms interact- 
ing with two intense pump fields which serve to modify 
the response of a weak probe field. The probe field is 
coupled simultaneously to both an electric and a mag- 
netic dipole transition. Included are local field correc- 
tions (LFCs) and chirality, which is explicitly shown to 
help produce a negative refractive index as well as tune 
the system between an amplifying (parametric) and ab- 
sorbing medium. 

In section [n] the theory needed to calculate the index 
of refraction in a four-level system is introduced which 
includes local field effects and magnetoelectric cross cou- 
pling. The numerical results for a dense atomic gas of 



neon are presented in section IIIII showing the possibility 
of obtaining a negative refractive index with low absorp- 
tion. The paper concludes in section ITVl 



II. THEORETICAL MODELING 

Fig. 1 shows the 4-level system under consideration 
that was originally introduced by Thommen et al. [11| . 
Although here, the magnetic component of the probe 
field (B) is explicitly included which couples to transi- 
tion |1) <r^ |2) with Rabi frequency fls = ^12 ■ B/h 
where mi2 is the corresponding magnetic dipole moment. 
The electric component of the probe field (E) couples to 
transition |3) -^ |4) with Rabi frequency D,e — ^34 ■ E/h, 
where ^34 is the electric dipole moment. Both transitions 
are nearly degenerate, allowing the probe field with fre- 
quency ijji2{= "^34) to simultaneously test both the per- 
meability and permittivity of the system. Additionally, 
two control lasers with frequencies lo^^ and CL'I'4 induce the 
magnetoelectric cross coupling that must be taken into 
consideration when solving for the material's response 
functions. The two control beams couple to transitions 
1) o |3) and |2) o |4) with Rabi frequencies J7i3 and 
SI24, respectively. The system's detunings are defined as 
5ij = uj^j — \{ei — ej)/h\ for which i,j = 1,2,3,4 and 
i ^ j. Ei denotes the energy of the corresponding atomic 
level. 

In the rotating wave approximation, the Hamiltonian 
describing the 4-level system is written as. 
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where aij = \i){j\ {i,j S 1 — J- 4) are level shift operators. 
Given the closed loop scheme presented, at least one of 
the Rabi frequencies must be allowed to be complex [1^ . 
Assuming a dense atomic medium, local field correc- 
tions (LFCs ) ar e included for which mean-field theory is 
utilized J2(]l423l |. The macroscopic electric fields in the 
Hamiltonian are therefore replaced by the microscopic, 
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FIG. 1. 4-level system under consideration. A probe field 
with frequency uji2 couples the electric (|3) <r^ |4)) and mag- 
netic (|1) -f^ |2)) transitions with Rabi frequencies Q,e and 
Hb, respectively. Two intense fields with Rabi frequencies 
fli3 and f224 couple the optical transitions (|1) ■(-> |3)) and 
(|2) •<->■ |4)) with the former being a two photon transition 
due to parity [ll[|. The corresponding wavelengths for Neon 
are A12 = 5.4 fim, A13 = 704 nm, and A24 = 352 nm with a rel- 
atively small energy difference of 0.16 MHz between the two 
transitions coupled to the probe field. The spontaneous decay 
rates are denoted by 721, 731, 742 and 743 while the energy 
levels are given in jL-coupling notation. 



local fields according to the Lorentz-Lorenz relation given 
by, 



£,LFc ^ ^ ^ :fip, 



(2) 



with P = ri(d), where 77 is the density of atoms and (d) is 
the expectation value of the electric dipolc moment under 
consideration. Similarly, in order to consider LFCs on 
the magnetic field one can write an analogous equation 
to that above: 
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with AI = ?/(?Ti), and (m) yielding the expectation value 
of the magnetic moment. The external magnetic field is 
taken to be iJ = i? = ±iE with the minus or plus sign 
referring to a^ or a^ circularly polarized light, respec- 
tively. 

A general equation for the polarization or magnetiza- 
tion of a medium should include cross-terms in order to 
account for the possibility of chirality, i.e., part of the 
polarization (magnetization) is induced by the magnetic 
(electric) field. In order to calculate the system's perme- 
ability and permittivity in terms of the electric, as-, and 



magnetic a^ polarizabilities, the theory previously intro- 
duced by Kastel et al. [la] is followed with the equations 
for the polarization and magnetization given by. 
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where q;eh and ckhe arc defined as the cross-coupling po- 
larizabilities due to chirality. 

The permittivity (e) and permeability (/i) are defined 
in relation to the susceptibilities such that e = 1 + 47rxe 
and ^ = 1 -I- 47rxm with Xe(Xm) being the electric (mag- 
netic) susceptibility of the system. Assuming the mate- 
rial's overall response to be isotropic, one then makes use 
of the general definitions for polarization and magnetiza- 
tion to calculate the susceptibilities and eventually the 
index of refraction. Hence, within isotropic chiral media 
we have the definitions, 
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(7) 

with ^EH = 47r2agjj and ^he = 47r^aHE [H- Using 
Eqs. (4) and (2) (Eqs. (5) and (3)), one is able to obtain 
the permittivity and permeability in terms of polarizabil- 
ities, yielding 
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aHae). 



Without the effects of chirality, the last term in paren- 
thesis of dH), dSl), and dTO]) would be left out. The first 
two terms in (jlOp that include only an (qe) are due to 
LFCs, and if chirality is not included then only the q^h 
(oe) term would be added to /x(e), however with chirality 
both terms appear in the denominator of the permittivity 
and permeability. The index of refraction, which is now 
dependent on the cross-coupling polarizabilities, is de- 
rived from the Helmholtz equation and given by [1^, [2J1 : 
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where n^ refers to cr* polarization. Note that a negative 
index of refraction does not necessarily require that both 
the electric and magnetic responses of the material be 
negative at the same time. This lessens the demand to 
induce a strong magnetic response which as previously 
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FIG. 2. Shown is the real part of the permeabiUty as a func- 
tion of probe field detuning for various coupling field ratios 
f^24/f2i3 = 2 (solid line), r224/i^i3 = 3 (dashed-dotted line), 
and ^24/^^13 = 10 (dotted line). The shaded area indicates 
the interval over which the real part of the permeability takes 
negative values. (See text for additional parameters.) 
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mentioned has been one of the biggest obstacles to over- 
come [HHIii. 

To characterize the medium's response, the ratio be- 
tween the real and imaginary parts of the index of re- 
fraction {n = n -\- in ) is used as a comparison between 
the phase velocity and the rate of absorption. Ideally, 
one aims for a negative real part of the index of refrac- 
tion along with a low imaginary part representing low 
absorption. The figure of merit is defined as [IJ) ll3 '■ 
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III. NUMERICAL RESULTS 

In order to solve for the dynamics of the system, the Li- 
ouville equation is numerically simulated, p = —i [H, p] + 
Tp, where Tp accounts any decay associated with spon- 
taneous emission. For the four-level system under inves- 
tigation the full density matrix solutions are given by, 



Pll = -2lm{{flB + Li2Pi2)p21 + (f^l3 + -^13Pl3)P3l} + 721P22 + 731P33 
P22 = -2lm{{ri*g + Ll2P2l)Pl2 + (f^24 + L24P24:)P'l2} " 721P22 +742^44 
P33 = -2Im{(fii3 + Li3P3i)/9i3 + {Qe + i34/534)/043} "731^33 + 743P44 
Pl2 = iSl2Pl2 + ii^B + Ll2Pl2){P22 - P\l) + i{^lZ + -^13Pl3)P32 - j(f^24 



Pl3 = i'5i3Pi3 +i(r2i3 

Pi4 = ihiPu + iiP-B - 

h-i = *<523P23 + «(^B - 
P24 = i'524P24 +«(r224 
/O34 = J(534P34 +i(f7_E - 
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il2P2l)Pl3 -i{^\-i 
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Li3Pi3)p21 + «(f^24 + i24P24)P43 ^ T-i^M 
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- i{n*g + Ll2P2l)pi4 - i{^E + L34P34)p23 - r24P24 
i(51l3 + Lj3P3i)pi4 - i{?t24 + L24P24)P32 - r34P34 
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where L12 = |"ii2| ^TTi-i/{3h) and Ly = |dyp47r7?/(3/i) 
with ij = {13,24,34}. These terms take into account 
the local field corrections which are applied to all of the 
fields within the Hamiltonian according to the definitions 
given in Eqs. 1^ and ^. The dipole moments are taken 
to be real, hence d^ = dji and TO12 = m2i. The decay 
of the off-diagonal elements of the density matrix, Fjj, 
are given in the radiative limit. In addition, we consider 
the Rabi frequencies of the coupling beams to be real 
while the Rabi frequencies of the probe field are taken to 
be complex with fls = fi^e*"^ and Qe = ^%e^'^. It has 
previously been shown that by adjusting the phase an- 
gle (j), one can minimize the imaginary part of the index 
of refraction, hence maximizing absorption or amplifica- 
tion of the probe field [1^, [13 ■ For our simulations, the 
special case of = tt is used. Also note that as seen in 
Eqs. (16)-(21), the LFCs yield the following effective de- 



tuning JJ^J'^ 
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Lij{pjj — Pa), known as the Lorentz 
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By numerically integrating the density matrix equa- 
tions, we are able calculate the material response for a 
wide range of parameters. The coherences between the 
states coupled to the probe field yield the polarization, 
P = ?7P43(i34, and magnetization. A/ = ?/P2i"ii2- Fol- 
lowing Rcf. [i5|, in order to calculate the polarizabilitics 
from Eqs. (4) and (5) we expand the coherences of the 
probe field transitions as power series in E^^"^ and H^^'^ . 
From the polarizabilitics the permittivity and permeabil- 
ity arc then obtained from Eqs. (8) and (9) and the index 
of refraction from Eq. (11). 

Figs. 2 and 3 display the real part of the perme- 
ability and pcrmitivity, respectively, for different ra- 
tios of fi24/f^i3 with (5i3 = —1.5 X lO^s^^, 77 = 1.0 X 
10-^^at/cm^,73i = 743 = 742 = 10'^s^\ and 721 = 
10^/a^ s~"^ where a = 1/137 is the fine structure con- 
stant. The largest Rabi frequency $^24 = 15 x 10^ s~^ 
is kept constant throughout the simulations correspond- 
ing to a laser intensity on the order of 10^ W/cm^. For 
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FIG. 3. The real part of the permittivity is plotted as a 
function of the probe detuning, S12, for various ratios be- 
tween coupling fields; r224/f^i3 = 2 (solid line), $124/^13 = 3 
(dashed-dotted line), and il24/r2i3 = 10 (dotted line). 



the case of Neon gas, this is well within the reported 
breakdown intensity of approximately lO^"*^ W/cm^ for a 
density on the order of 10'^^ atoms/cm^ if one is using a 
pulse of 7ps[23- The strong effects of chirality begin to 
be seen as the ratio between coupling fields is increased 
and we notice the appearance of multiple peaks in the 
permeability and permitivity. This is due to the depen- 
dence of the cross-polarizabilities not only on the real 
parts but also the imaginary parts of the coherences of 
the probe field transitions. At very large dctunings, the 
material's response functions (not shown on Figs. 2 and 
3) converge towards the expected values for a field trav- 
elling through free space. 

Fig. 4 shows the calculated results for the index of 
refraction which becomes negative for an interval of ^ 
50 X 10^ s~^ (shaded area) over the detuning of the probe 
field with a figure of merit of roughly 30 corresponding 
to a coupling ratio of ^224/^13 = 10 (see Fig. 4(b)). Note 
that this range of probe field detuning exhibits a nega- 
tive index of refraction larger than the range in which the 
permitivity (shaded area in Fig. 2) is negative. The real 
part of the two chiral terms ^eh and ^he (inset) plays 
a crucial role in decreasing absorption as the difference 
between the two directly affects the imaginary part of the 
index of refraction (see Eq. (fTT|) ). The two peaks corre- 
sponding to the Re(^H£;) help produce the two troughs 
in the curve representing n and yield the highest val- 
ues for the figure of merit, thus providing the points for 
which we find the minimum amount of absorption. By 
adjusting the ratio between coupling fields ^24 and ^13, 
we have ourselves a useful tool with which we can fur- 
ther increase the influence of the second term outside the 
square root of Eq. (11) in order to manipulate the index 
of refraction. 

More recently, research has focused on producing 
a negative refractive index under conditions of trans- 
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FIG. 4. Shown are (a) the real (solid line) and imaginary 
(dashed line) parts of the index of refraction, and (b) the 
values for the figure of merit as a function of probe detuning. 
A negative index of refraction is shown to roughly occur for a 
detuning between 14.3xl0*s~^ and 15x10* s~^ (shaded area) 
with a FoMr;30. Results presented are for a coupling ratio 
of ^24/^^13 = 10. The inset in (a) shows the real part of the 
chiral terms which contribute to the imaginary part of the 
index of refraction. 



parcncy or gain [a lul, l26|, l27|, l29|, |30|. Utilizing the 
method mentioned above, the ratio between the Rabi fre- 
quencies of the two coupling beams is adjusted in order 
to increase chirality and produce the desired transparent 
or amplifying effect. Plotted in Fig. 5, both the real (a) 
and the imaginary (b) part of the index of refraction is 
shown as a function of probe field detuning for different 
values of this ratio. In all the cases, one obtains a neg- 
ative index of refraction along with a decrease in the n 
as the ratio between the two coupling fields is increased. 
For the value r224/^i3 = 20 the imaginary part of the 
index of refraction becomes negative for a range of ap- 
proximately 20 X 10® s~^ over the probe field detuning 
signifying probe gain. This gain can be directly related 
to a chiral enhancement of the probe field's electric tran- 
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FIG. 5. Shown is (a) the real and (b) the imaginary part of the 
index of refraction as a function of the probe detuning, (5i2, for 
coupling field ratios of Vl^aJ^iz = 10 (solid line), 024/^13 ~ 
15 (dashed-dotted line), and f224/f2i3 = 20 (dotted line). 



FIG. 6. Shown is (a) the real and (b) the imaginary part of 
the index of refraction as a function of the probe detuning, 
(5i2, for different values of the atomic density (r\). All other 
parameters are identical to the case for neon gas presented in 
Fig. 4. 



sition [28|. 

The role of LFCs in the response of the medium 
has also been herein investigated and previously demon- 
strated to play a significant role in inducing transparency 
[20|. The effect of local field corrections can easily be 
regulated by changing the atomic density, r\. As shown 
in Fig. 6, as the number of atoms, and therefore the 
strength of the local field corrections, is increased we 
see a corresponding decrease in absorption as well as 
the expected Lorentz shift. In particular, for the value 
77 = 1.5 X lO^^at/cm'' we obtain a range of detunings for 
which the medium begins to amplify. 



IV. CONCLUSIONS 

In this paper, a method for producing a negative in- 
dex of refraction for a probe field within the infrared fre- 
quency of light is presented. A specific, practical model 
was demonstrated representing a dense Neon gas which 
has two near degenerate electric and magnetic transi- 
tions. Taken into consideration was the presence of cross 



polarizability terms in the material's response functions 
due to a strong chirality which can be utilized to ma- 
nipulate the refractive index. It has been shown that by 
adjusting the ratio of two intense pump fields, which cou- 
ple a nearly degenerate electric and magnetic transition, 
one can manipulate the effects of chirality in order to in- 
duce low absorption or gain in the system. In addition, 
local field corrections were added to the calculation where 
it was shown that the atomic density could also be used 
as a tool to further reduce absorption. The use of gases 
has certain potential benefits compared to metamateri- 
als, for instance, eliminating manufacturing contraints 
as well as reducing experimental difficulties such as the 
alignment of laser and material polarizations. However, 
it is important to note that these two methods that have 
been employed for producing a negative refractive index 
with low absorption can be used not only in a dense 
atomic gas, but within similar metamaterial structures 
exhibitingthe inherent traits of chirality and requiring 
LFC'sQ, |31|. The desired result, a consistent source of 
negatively refracted light, is crucial for a multitude of 



applications, particularly lenses with perfect resolution, 
subwavelength enhancement, as well as the generation 
transparent materials 0, [S^. Although this paper has 
considered an isotropic, infrared source there are no re- 
strictions to use a material of this nature and indeed 
further research is being undertaken to potentially show 
how a non-isotropic light source could be used as an ad- 
ditional tool in adjusting the Helmholtz equation and 
allowing increased control over the refractive index. 
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